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Abstract 

The absolute muon flux between 20 GeV and 3000 GeV is 
measured with the L3 magnetic muon spectrometer for zenith 
angles ranging from 0° to 58°. Due to the large exposure of 
about 150 m^ sr d, and the excellent momentum resolution of the 
L3 muon chambers, a precision of 2.3 % at 150 GeV in the vertical 
direction is achieved. 

The ratio of positive to negative muons is studied between 
20 GeV and 500 GeV, and the average vertical muon charge ratio 
is found to be 1.285 ± 0.003 (stat.) ± 0.019 (syst.). 

The L3+C group dedicates this publication to the late Bianca Monteleoni. 



Submitted to Phys. Lett. B 



Introduction 



Atmospheric muons are among the final products of cosmic ray induced air- 
showers. The absolute muon flux and its momentum dependence are mainly 
determined by the flux of nucleons entering the atmosphere and the inclusive 
meson production cross sections in high-energy hadronic interactions. The 
ratio of the fluxes of positive to negative muons, denoted as charge ratio in 
the following, reflects the proton to neutron ratio at the top of the atmo- 
sphere, folded with the production and decay spectra of charged pions and 
kaons. While the knowledge of the primary cosmic ray spectrum below a 
few 100 GeV has improved considerably in the recent past [1], large uncer- 
tainties still exist in the primary energy range between 0.1 TeV and 500 TeV 
responsible for the production of secondaries with momenta in the range un- 
der study here. Moreover, the details of high energy hadronic interactions 
still lack theoretical understanding and there is little experimental data in 
the relevant energy and phase space regions [2]. Therefore the ground- level 
muon flux and charge ratio are widely used to tune or verify the parameters 
of atmospheric cascade calculations [3-6]. Currently these calculations are 
of great interest, as they predict the absolute atmospheric neutrino fluxes [7] 
which are needed to interpret the observed muon neutrino flux deficit [8-10] 
and to evaluate the backgrounds for neutrino astronomy. 

The muon fiux and charge ratio have been extensively studied with dif- 
ferent experimental techniques [11]. However, results show discrepancies of 
about 10—20% with respect to each other, which exceed the uncertainties 
assigned to the individual measurements and thus indicate the presence of 
systematic effects not accounted for. 

Here a new measurement of the atmospheric muon fiux is presented using 
the precise muon spectrometer of the L3 detector located at the LEP collider 
at CERN, near Geneva, Switzerland. Special attention is given to the precise 
determination of all relevant detector and environmental parameters needed 
to convert the raw-data distributions into an absolute surface level fiux. The 
large statistics available permits extensive studies of the residual systematic 
uncertainties. 

Experimental setup 

The momentum distribution of atmospheric muons is measured with the 
upgraded L3 setup of the L3 detector [12] known as L3-I-C [13]. The parts 
of the detector used in this analysis are sketched in Figure 1. After passing 
through the stratified rock overburden, called "molasse", the arrival time to 
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of a muon is measured with a resolution of 1.7 ns by a 202 scintillator 
array placed on top of the L3 detector. The array is composed of 34 modules, 
each read out by two photomultipliers in coincidence to reduce noise. Inside 
a volume of about 1000 m'^ with a magnetic field of 0.5 T, the coordinates and 
slopes of a muon track arc measured in up to six drift chambers in the bending 
plane and up to eight times in the non-bending plane. These chambers are 
arranged concentrically around the LEP beam line in two groups of eight 
octants, each containing three layers of drift cells. By subtracting the 
time from the arrival times of the drift electrons at the sense wires, a track 
position in each chamber can be reconstructed with a precision of about 
60 /im in the bending plane and 1 mm in the non-bending plane. 

Only three points are needed to determine the radius of the track in the 
magnetic field, therefore the momentum of a muon traversing two octants 
can be measured twice. This redundancy is used to evaluate the detector ef- 
ficiencies and the resolution of the apparatus. The best resolution is obtained 
when fitting the six points together over the full track length of 11 meters. 
The multiple scattering and energy loss inside the L3 inner detectors, as well 
as the effect of the inhomogcncous magnetic field are taken into account using 
the procedure proposed in Reference [14]. 

Equipped with a trigger and data acquisition system independent of the 
normal L3 data-taking, L3-I-C recorded 1.2 x lO^'^ atmospheric muon triggers 
during its operation in the years 1999 and 2000. 

The L3+C experiment was located 450 m above sea level at a longitude 
of 6.02° E and a latitude of 46.25° N. 

For vertically incident muons, the mean energy loss in the molasse {X— 
6854 gcm~^) and the magnet (X=1227 gcm~^) is 19 GeV at low momenta 
and reaches 57 GeV at 1 TeV. 

Analysis 

Detector and molasse simulation 

The geometrical acceptance of the L3-I-C detector and the stochastic energy 
loss in the molasse overburden are evaluated using the following simulation 
procedure: Monte Carlo events are generated on the surface using a parame- 
terization of the zenith angle and momentum dependence of the muon spec- 
trum as obtained with the CORSIKA [15] program. These simulated muons 
are then tracked through a GEANT [16, 17] model of the L3-I-C environment 
which includes the molasse, access shafts and the concrete structures around 
the cavern which hosts the apparatus. Finally, the detector response is sim- 
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ulated with a detailed GEANT description of the L3 detector. The generated 
detector signals are reconstructed with the same program used for the data. 
In total 1.7 X 10^ reconstructed Monte Carlo events are used in this analysis. 

Event selection 

The data analysis is restricted to events with three position measurements 
in at least one octant, a scintillator hit and good running conditions during 
data-taking. A total of 1.2 x 10^ reconstructed muon tracks are retained. 

The shielding of the 30 m of molasse overburden absorbs most of the 
charged air-shower particles other than muons. The number of muons pro- 
duced in e'^e" collisions by LEP is negligible compared to the flux of at- 
mospheric muons. Therefore no background rejection is needed. The data 
selection focuses on two topics. Firstly, fiducial volume cuts are defined to 
assure a good description of the data by the simulation. Secondly, selection 
cuts are imposed on the track quality to enhance the momentum and angular 
resolution. These selection criteria are: 

• The muon track positions must be measured in six layers in the bending 
plane. 

• The momentum resolution, calculated from the quality of the track 
position measurements, should not exceed its nominal value by more 
than 50%. 

• At least four position measurements (two in each octant) should be 
present in the non-bending plane. 

• The of a fit of the tracks to a circle within an octant must satisfy 
xVndf < 4. 

• The difference between the two photomultiplier time measurements 
from the same scintillator module must be below 8 ns. 

After these cuts, 2 x 10^ data events remain for the muon spectrum analysis. 
The selection efficiencies depend on the charge, momentum and direction of 
the muon. For muons with momenta above 80 GeV, the average efficiency is 
7.6 % for the fiducial volume cut and 33.3 % for the quality selection. 
4 % of the raw events are multi-muon events. Each individual muon is 
counted as an input to the spectrum data. 



4 



Momentum resolution 

The single-octant resolution is inferred directly from the data by comparing 
the two independent curvature measurements of muons traversing two oc- 
tants. An example of the curvature difference is shown in Figure 2. These 
measurements are used to tune the detector simulation, from which the reso- 
lution of the full fit is determined. The relative momentum resolution Ap/p 
as a function of momentum at the detector- level is shown in Figure 3(a). The 
maximum detectable momentum of the spectrometer, defined as the momen- 
tum at which Ap/p reaches unity, is 0.78 TeV for muons measured in only 
one octant and about 5 TeV for muons measured in two octants. 

Detector efficiencies 

The efficiency of each subdetcctor is stTidicd by exploiting redundancies in 
the measurement process. For about 50 % of the tracks, the muon arrival 
time is also deduced from the muon chambers. These tracks arc used to 
determine the scintillator efficiencies as a function of time and position on 
the array. A mean efficiency of 95.6% is found at the start of data-taking 
decreasing continuously to 94.5 % towards the end of 2000. The possibility of 
reconstructing a muon within a single octant is used to scan the drift-layer 
performance of the facing octant. On average, a fraction of 10.5% of the 
drift cells are found to have an efficiency lower than 80 %. These regions 
are excluded in both the data and Monte Carlo reconstruction. Under these 
conditions the trigger efficiency is determined from redundant trigger classes 
to be 99.85 % on average. 

During data-taking, the total effective running time was continuously 
measured with a 10 MHz live-time counter, which is disabled whenever the 
trigger system is not ready to accept new data. In addition, each second an 
external trigger signal was sent to the L3-I-C trigger system. The number 
of these external triggers on tape compared to the total number of running 
seconds gives another estimate of the effective running time and agrees with 
the value from the live-time counter within 0.02 %. 

Selection efficiency 

Using the possibility to measure a muon independently in two detector parts, 
the selection efficiencies are determined in the following way: the detector 
is subdivided into two hemispheres, i and j, and the conditional hemisphere 
selection probabilities are measured for data and Monte Carlo separately 
as a function of the muon charge g, momentum p and zenith angle 9. In the 



5 



absence of correlated inefficiencies, tfie total selection efficiency for accepting 
a track in the two hemispheres is given by the product £i x £2- The ratio 



r = 



(1) 



is used to correct the differences between data and Monte Carlo. Depending 
on the zenith angle range and the data-taking year, r varies from 0.84 to 
0.90. A large fraction of this correction factor originates from a defect in the 
TDCs used to read out the muon chambers, giving rise to an 8 % inefficiency 
for the full track selection. 

Surface spectrum 

The relation between the momentum distribution measured in L3-I-C and the 
muon surface spectrum is given by 



Here n is the vector of events rii with measured momenta between [qp-, Q'Pj+i]. 
The effective live-time is given by r and R denotes the migration matrix, i. e. 
the conditional probability of measuring a momentum qp^ given a surface 
momentum qpj. A is the diagonal matrix of geometrical acceptances as a 
function of the surface momentum and E is the diagonal matrix of detector 
efficiencies as a function of momentum at the detector-level. The vector 
m contains the true surface spectrum integrated over a surface momentum 
bin. The complete detector matrix, D = E • R • A, is evaluated from the 
measured detector efficiencies and the detector simulation as follows: 



where Sue is the surface area used in the Monte Carlo generator, AQ the solid 
angle of the zenith bin under study, £j includes the scintillator and trigger 
efficiencies and is the selection efficiency correction discussed above, n^^' 
denotes the number of selected Monte Carlo events found within a detector- 
level momentum bin i, which were generated within the momentum bin j at 
the surface, and N^^^ is the total number of Monte Carlo events generated 
within this surface momentum bin. 

The effective acceptance of this analysis is calculated by summing over 
the columns of the detector matrix, which yields the geometrical factor for a 
muon being registered in the detector and fulfiUing the selection cuts. The 



n = TERAm. 



(2) 




(3) 
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product of the effective acceptance and tfie live-time gives tfie total expo- 
sure, shown in Figure 3(b) for positive and negative muons as a function 
of surface momentum. It rapidly decreases at low energies due to the mo- 
mentum cut-off caused by the molasse overburden. Below 200 GeV, positive 
and negative muons have different acceptances, because the magnetic field 
bends their tracks in opposite directions and correspondingly into different 
detector regions. At large momenta the acceptance decreases with the per- 
formance of the full detector fit and a more difficult reconstruction caused 
by the increasing production of delta rays. 

The measurement Equation (2) is solved using the least squares method, 
by minimizing 

^ V <7,(m)2 ' 

where cr, contains the statistical errors of the data and the detector matrix: 



(7,; = 




J^Vmml (5) 



In the first step, the statistical Monte Carlo variances are set to zero, 

such that Equation (4) becomes hnear with respect to the surface spectrum m 
and its solution is 

m= ^ (D^WD)"^D^Wn (6) 

with covariance matrix 

V[m] = (D^^WD)"^ . (7) 

Here W denotes the diagonal weight matrix containing the statistical errors 
of the data and the Monte Carlo, Wa = l/crf. 

The minimization of Equation (4) is then repeated using the solution m 
of the previous iteration for the calculation of the errors in Equation (5). 
This process is repeated until the maximum relative difference to the result 
of the previous iteration is below 10~^. Typically four iterations are needed. 



Systematic uncertainties 

Normalization uncertainties 

Uncertainties on the live-time and the trigger and scintillator efficiencies give 
rise to a normalization uncertainty of 0.7%. 
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The uncertainty of the detector acceptance is assessed in three studies: 
First, the results obtained for statistically independent data subsamples, as 
for instance the data collected in 1999 and 2000 or in different detector parts, 
are compared. Second, the muon flux and charge ratio are measured as a 
function of the azimuthal angle. At large momenta, geomagnetic effects and 
the variation of the molasse overburden are not important, and therefore a 
flat distribution is expected. Finally, the stability of the measured flux and 
charge ratio with respect to a variation of the selection criteria is investigated. 

From these studies, additional normalization uncertainties in the absolute 
muon flux are derived. These range from 1.7% to 3.7% depending on the 
zenith angle. For the charge ratio normalization uncertainties between 1.0% 
and 2.3% are estimated. 

Figure 4 shows an example of the stability of the measured muon flux 
within the two data-taking years for zenith angles between and 32 de- 
grees. Since the muon production is known to change with atmospheric 
conditions [18, 19], the observed variation of the muon flux is compared to an 
air-shower simulation with the TARGET [20] program using atmospheric den- 
sity proflles measured in balloon flights close to the experiment [21]. Good 
overall agreement between data and Monte Carlo is observed. However, the 
full comparison to the rates in 26 weeks and 14 momentum bins yields a 
X^/ndf of 526/364. The assumption that this large value is caused by de- 
tector inefficiencies not accounted for, leads to an additional normalization 
uncertainty of 0.3 %, which is well within the above estimated uncertainties. 



Momentum scale uncertainties 

Due to the steepness of the muon spectrum, even small uncertainties in the 
absohite momentum scale can introduce a considerable bias in the muon flux 
measurement. 

The uncertainty on the L3 magnetic field strength introduces a momen- 
tum scale bias of less than 0.4% [22]. 

Furthermore, the momentum measurement is subject to uncertainties of 
the detector alignment. A systematic shift of the chamber positions may 
introduce a constant offset C . The measurement of the curvature, g/p, and 
the alignment related momentum scale uncertainty, 5ai is given by 

AC 

and depends on the muon charge. Within one octant, the alignment is mea- 
sured by an optical alignment system [23] with a precision corresponding to 
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0.19 TeV-^ [24]. The relative alignment of the muon chamber octants, rel- 
evant for this analysis, is determined from the data itself with a precision 
between 0.075 and 0.152 TeV~^ [25], depending on the zenith angle. 

The Tinccrtainty due the molasse overburden affects the conversion of 
the measured flux at the detector to the surface. The results of two survey 
drillings at different locations close to L3-I-C provide an absolute measure- 
ment of the L3+C matter overburden. The influence of molasse inhomo- 
geneities and of surface installations not included in the L3+C simulation 
is estimated by studying the variance of the muon flux as a function of the 
azimuthal angle near the momentum threshold. This leads to an uncertainty 
of the average rock density of 2%, which is equivalent to an energy loss 
uncertainty of 0.4 GeV in the vertical direction. 

Good agreement between the muon energy-loss calculation used here [16, 
17] and other approaches [26-28] is found. The residual differences corre- 
spond to a momentum scale uncertainty below or less than 0.3 % in the 
vertical direction. 

The relative momentum scale uncertainties for vertically incident muons 

are displayed in Figure 5(a). At low energies the molasse uncertainty con- 
tributes the most, whereas above 100 GeV the alignment uncertainties dom- 
inate. 

Detector matrix uncertainty 

The limited Monte Carlo statistics affects the precision of the detector matrix 
D. Below 200 GeV, it dominates the total statistical uncertainty in the de- 
nominator of Equation (4), contributing about 0.5 % to the total uncertainty 
per zenith angle bin. 

In order to estimate the influence of the uncertainty of the momentum 
resohition on the measured muon flux, the minimization of Equation (4) is 
repeated with different detector matrices, for which the momentum resolu- 
tion is altered by ±8%. This corresponds to the estimated uncertainty of 
its Monte Carlo prediction. As expected, no differences are found at low 
momenta. Above 200 GeV, the observed relative flux change A$ is well 
described by 

A$ = c ■ {p-0.2 TcV), (9) 

with c = 0.03 TcV^^. The observed difference between the high-energy muon 
flux measured in different detector regions leads to a somewhat larger value 
of c = 0.06 TeY-\ 
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Total uncertainty 

The total uncertainties of the muon flux and charge ratio are obtained by 
adding the individual contributions in quadrature. The different sources of 
the vertical uncertainties are shown in Figures 5(b) and (c). The muon flux 
uncertainty is dominated by the uncertainty of the molasse overburden at low 
momenta and by the alignment and resolution uncertainty at high momenta. 
The minimal uncertainty is 2.3% at 150 GeV in the vertical direction. The 
vertical charge ratio uncertainty is below 2% up to momenta of 100 GeV. 
Above this momentum, it rises rapidly with the alignment uncertainties. 

These uncertainties are fully correlated between different momenta for a 
given zenith angle bin. As approximately the same detector parts arc used to 
measure the muons in neighboring zenith angles, the systematic uncertainties 
are also correlated with respect to the zenith angle. The estimated correlation 
coefficients are hsted in Table 1. 

Z-events 

The understanding of the detector is validated by analyzing the muons pro- 
duced at LEP via the process 

e"'"e~ — >• Z — >• , 

recorded during the LEP calibration runs at a mean centre-of-mass energy 
of 91.27 GeV. The selection criteria include the requirement of a muon track 
close to the collision point and an event-time in coincidence with the LEP 
beam crossing time. The number of selected muons with a momentum above 
60 % of the beam energy is converted to an absolute cross section resulting 
in 

(7j;3+c = 1.447 ± 0.071 (stat.) ± 0.021 (syst.) nb . (10) 

Here the quoted systematic uncertainty includes only sources which are not 
relevant to the muon spectrum measurement, such as the luminosity. 
Using the LEP precision measurements [29], the Standard Model expectation 
of o-fj^+fj^^ is calculated [30] to be 

a^f^. = 1.4840 ± 0.0013 nb , (11) 

which is in excellent agreement to the value measured here. 
Thus this study verifies the L3+C acceptance calculation and a normaliza- 
tion uncertainty of <5.2% at 68% C.L. can be stated. Although this number 
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is larger than the estimated systematic uncertainty of the muon flux normal- 
ization, it provides an absolute systematic cross-check qualitatively different 
from the relative studies described above. 

The momentum distribution of the selected events, displayed in Figure 6, 

shows good agreement between the data and the simulation. From the peak 

position of the data, an absolute momentum scale uncertainty of < 370 MeV 

and a single octant alignment uncertainty of < 0.1 TeV~^ is derived. 

As can be seen in figure Figure 3(a), the momentum resolution derived from 

the Z muon sample agrees well with the one measured with atmospheric 

muons. 



Results 

The muon fluxes, conventionally multiplied by the third power of the 

momentum, and the charge ratios, i?, are listed for each zenith angle bin in 
Tables 2—9 with their statistical and systematic uncertainties. The statistical 
correlation coefficients, p, between neighboring momentum bins, as derived 
from Equation (7), are also given. Due to the limited detector resolution 
these correlations are inevitable. However, the momentum binning is chosen 
such that only neighboring bins have a significant correlation. 

The average momenta, (p), within a momentum range [j5i,p2] are calcu- 
lated [34] by fitting the phenomenological muon fiux function from Refer- 
ence [11] to our data and solving 

^((p)) = / ^{p)dv. (12) 

It should be noted that the fiuxes are neither corrected for the altitude 
of L3-I-C nor for the atmospheric profile to avoid additional theoretical un- 
certainties. Instead, we quote the average atmospheric mass overburden X 
above L3-I-C, which was continuously measured with balloon flights from 
close to the experiment to altitudes of over 30 km [21]. The parameteriza- 
tion of Reference [35] is used to describe the mass profile X in g cm~^ as a 
function of the altitude h in km above sea level: 

^ ^ fA(/i;,-/i)("+i), h<ll 
X(h) = 1 ^ ' - 13 

[Be'^, h>ll 

A fit to the live-time weighted balloon data yields A — 8.078 x 10~^, B — 
1332, h = 39.17, ho = 6.370 and a = 3.461. 

The measured muon fiuxes at the L3-I-C altitude are shown in Figure 7 for 
each zenith angle bin. As no previous continuous zenith angle measurements 
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exist in the large energy range examined here, only the vertical flux can be 
compared to other experiments, as shown in Figure 8. Only measurements 
providing an absolute normalization [36-44] are taken into account. The 
data are extrapolated to sea level using the muon flux predictions of the 
TARGET [20] program. 

The comparison to low energy experiments [38-43] gives a good overall 
agreement with this analysis above about 40 GeV. At lower momenta, a 
systematic slope difi^erence seems to be present, which corresponds to about 
three standard deviations of the systematic molasse uncertainty estimated 
above. 

Only three previous experiments measured a normalized spectrum at high 
energies. The shape of the Kiel measurements [36] agree with this result over 
the full momentum range, but a lower flux normalization is determined by 
L3+C. 

The data obtained with the MARS apparatus [37] significantly disagree with 
this result, both in shape and normalization. 

Above momenta of 50 GeV, the recent muon flux measurement from BESS- 
TeV [44] is in very good agreement with this result. 

The measured charge ratios at the L3+C altitude are shown in Figure 9 
for each zenith angle bin up to momenta of 500 GeV. In the considered 
momentum range, the charge ratio is independent of the momentum within 
the experimental uncertainties. The mean value in the vertical direction is 
found to be 1.285 ± 0.003 (stat.) ± 0.019 (sys.) with a xVndf =9.5/11. This 
is in good agreement with the average of all previous measurements, 1.270 
± 0.003 (stat.) ± 0.015 (sys.) [11]. It is worth noting, that the precision of 
the data of a single L3+C zenith angle bin is comparable to the combined 
uncertainty of all data collected in the past. 
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COS 6* 


0.525- 
0.600 


0.600- 
0.675 


0.675- 
0.750 


0.750- 
0.825 


0.825- 
0.900 


0.900- 
0.938 


0.938- 
0.975 


0.975- 
1.000 


0.525-0.600 


1.00 


0.91 


0.99 


0.94 


0.93 


0.00 


0.00 


0.00 


0.600-0.675 


0.91 


1.00 


0.96 


0.76 


0.88 


0.00 


0.00 


0.00 


0.675-0.750 


0.99 


0.96 


1.00 


0.91 


0.95 


0.00 


0.00 


0.00 


0.750-0.825 


0.94 


0.76 


0.91 


1.00 


0.93 


0.00 


0.00 


0.00 


0.825-0.900 


0.93 


0.88 


0.95 


0.93 


1.00 


0.01 


0.01 


0.01 


0.900-0.938 


0.00 


0.00 


0.00 


0.00 


0.01 


1.00 


1.00 


1.00 


0.938-0.975 


0.00 


0.00 


0.00 


0.00 


0.01 


1.00 


1.00 


1.00 


0.975-1.000 


0.00 


0.00 


0.00 


0.00 


0.01 


1.00 


1.00 


1.00 



Table 1: Correlation coefficients of the detector-related systematic uncertain- 
ties between different zenith angle bins from 0° to 58° 



19 



momentum 


ip) 


* ■ ip)' 




P4> 




R 


A stat 


PR 


^R 


interval 


[GcV] 


r Q^,v2 ] 


[%] 




l%] 




l%] 




l%] 


[vjtiv J 




















20.0-27.0 


23.18 


0.217 


0.4 


-0.24 


4.4 


1.274 


0.7 


-0.24 


1.5 


















27.0-34.5 


30.47 


0.231 


0.3 


-0.22 


3.7 


1.284 


0.6 


-0.22 


1.5 


















34.5-42.0 


38.02 


0.244 


0.4 


-0.24 


3.3 


1.295 


0.8 


-0.24 


1.5 


















42.0-50.0 


45.78 


0.252 


0.5 


-0.26 


3.0 


1.269 


0.9 


-0.26 


1.5 


















50.0-58.5 


54.04 


0.257 


0.5 


-0.31 


2.8 


1.286 


1.0 


-0.31 


1.5 


















58.5-68.5 


63.25 


0.261 


0.5 


-0.31 


2.6 


1.298 


1.0 


-0.31 


1.5 


















68.5-81.5 


74.63 


0.266 


0.5 


-0.30 


2.5 


1.273 


1.1 


-0.29 


1.5 


















81.5-100 


90.13 


0.262 


0.6 


-0.26 


2.4 


1.293 


1.2 


-0.26 


1.5 


















100-132 


114.5 


0.264 


0.5 


-0.22 


2.3 


1.286 


1.0 


-0.21 


1.6 


















132-200 


161.3 


0.250 


0.5 


-0.22 


2.2 


1.287 


1.1 


-0.22 


2.1 


















200-300 


243.0 


0.233 


0.8 


-0.27 


2.2 


1.327 


1.7 


-0.27 


3.7 


















300-500 


381.9 


0.203 


1.2 


-0.28 
-0.30 


2.4 


1.276 


2.6 




7.2 




















500-1000 


687.2 


0.151 


2.3 


3.5 










1000-3000 


1599 


0.087 


6.0 




8.4 











Table 2: Muon flux, $, multiplied with the third power of the momentum, 
and charge ratio for 0.975 < cos^ < 1.000. The statistical, A^*^*, and the 
systematical, A^^^^, uncertainties are given. p$ and pr are the statistical 
correlation coefficients between neighboring momentum bins, as derived from 
Equation (7). 



momentum 


(p> 


* ■ (P>^ 


A stat 


P4 




R 


A stat 
^R 


PR 




interval 


[GeV] 


r GeV^ 1 
[ cm-^s sr J 


[%] 




[%] 




[%] 




[%] 


[GeV] 




















27.0-34.5 


30.47 


0.237 


0.3 


-0.23 


4.5 


1.265 


0.6 


-0.22 


2.3 


















34.5-42.0 


38.02 


0.254 


0.4 


-0.24 


4.1 


1.261 


0.8 


-0.24 


2.3 


















42.0-50.0 


45.78 


0.259 


0.5 


-0.26 


3.8 


1.265 


1.0 


-0.26 


2.3 


















50.0-58.5 


54.04 


0.264 


0.5 


-0.30 


3.7 


1.297 


1.1 


-0.30 


2.3 
















58.5-68.5 


63.25 


0.266 


0.6 


-0.31 


3.5 


1.250 


1.1 


-0.31 


2.3 


















68.5-81.5 


74.63 


0.270 


0.6 


-0.29 


3.4 


1.319 


1.2 


-0.29 


2.3 


















81.5-100 


90.13 


0.271 


0.6 


-0.25 


3.3 


1.259 


1.3 


-0.25 


2.3 


















100-132 


114.5 


0.270 


0.5 


-0.21 


3.2 


1.296 


1.1 


-0.21 


2.5 


















132-200 


161.3 


0.262 


0.6 


-0.21 


3.1 


1.281 


1.2 


-0.21 


3.0 


















200-300 


243.0 


0.241 


0.9 


-0.27 


3.1 


1.273 


1.8 


-0.27 


4.8 


















300-500 


382.1 


0.219 


1.3 


-0.27 


3.2 


1.389 


2.7 




8.3 




















500-1000 


687.8 


0.165 


2.4 


-0.30 


4.2 




























1000-3000 


1604 


0.093 


6.5 




9.2 











Table 3: Muon flux and charge ratio for 0.938 < cos6l < 0.975. 
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momentum 




* ■ {pf 




P4 




R, 


A S t tit 

R 


PR 


^R 


mtcrval 


[GcV] 


r 1 n 


[%] 






[%] 




[%] 




[%] 


[GeV] 




[ cm-^s sr J 


















27.0-34.5 


30.47 


0.245 


0.4 


-0 


24 


4.3 


1.265 


0.7 


-0.24 


2.3 


















34.5-42.0 


38.02 


0.257 


0.4 


-0 


26 


3.9 


1.296 


0.9 


-0.26 


2.3 


















42.0-50.0 


45.78 


0.265 


0.5 


-0 


28 


3.6 


1.302 


1.1 


-0.28 


2.3 


















50.0-58.5 


54.04 


0.273 


0.5 


-0 


31 


3.4 


1.267 


1.1 


-0.31 


2.3 


















58.5-68.5 


63.25 


0.279 


0.5 


-0 


31 


3.2 


1.287 


1.1 


-0.30 


2.3 


















68.5-81.5 


74.64 


0.282 


0.5 


-0 


27 


3.1 


1.272 


1.1 


-0.27 


2.3 


















81.5-100 


90.14 


0.283 


0.5 


-0 


23 


3.0 


1.288 


1.1 


-0.23 


2.4 


















100-132 


114.5 


0.284 


0.5 


-0 


19 


2.9 


1.303 


0.9 


-0.19 


2.6 


















132-200 


161.3 


0.274 


0.5 


-0 


19 


2.8 


1.325 


1.0 


-0.19 


3.8 


















200-300 


243.0 


0.255 


0.7 


-0 


23 


2.8 


1.300 


1.4 


-0.23 


6.8 


















300-500 


382.0 


0.220 


1.1 


-0 


25 


2.9 


1.437 


2.2 




12.0 




















500-1000 


687.3 


0.172 


1.9 


-0 


29 


4.2 




























1000-3000 


1599 


0.091 


5.7 






11.3 











Table 4: Muon flux and charge ratio for 0.900 < cos 6' < 0.938. 



momentum 


(p> 


* ■ (pf 




P4 




R 


A stat 
^R 


PR 




interval 


[GeV] 


r GeV^ 1 
1 cm-^ s sr J 


[%] 






[%] 




[%] 




[%] 


[GeV] 






















34.5-42.0 


38.02 


0.263 


0.4 


-0 


28 


3.9 


1.248 


0.8 


-0.28 


1.3 


















42.0-50.0 


45.78 


0.269 


0.5 


-0 


29 


3.5 


1.288 


1.1 


-0.29 


1.3 


















50.0-58.5 


54.04 


0.274 


0.5 


-0 


32 


3.3 


1.280 


1.1 


-0.32 


1.3 


















58.5-68.5 


63.25 


0.284 


0.6 


-0 


31 


3.1 


1.285 


1.2 


-0.31 


1.3 


















68.5-81.5 


74.64 


0.280 


0.6 


-0 


28 


2.9 


1.274 


1.2 


-0.28 


1.3 


















81.5-100 


90.14 


0.288 


0.6 


-0 


24 


2.8 


1.302 


1.3 


-0.24 


1.3 


















100-132 


114.5 


0.290 


0.5 


-0 


20 


2.6 


1.280 


1.1 


-0.20 


1.4 


















132-200 


161.3 


0.279 


0.6 


-0 


20 


2.6 


1.313 


1.1 


-0.20 


2.0 


















200-300 


243.1 


0.264 


0.8 


-0 


27 


2.5 


1.301 


1.7 


-0.27 


4.7 


















300-500 


382.1 


0.232 


1.3 


-0 


30 


2.6 


1.318 


2.6 




11.3 




















500-1000 


686.8 


0.167 


2.5 


-0 


37 


4.2 




























1000-3000 


1587 


0.087 


7.9 






10.4 











Table 5: Muon flux and charge ratio for 0.825 < cos 6* < 0.900. 
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momentum 




^ \t'i 


A Stat 


p* 


A sy.st 


R 


A stat 


PR 


A sy.st 


interval 


[GcV] 


r 9 1 

GoV^ 


[%] 




[%] 




[%] 




[%] 


[GeV] 




















34.5-42.0 


38.02 


0.255 


0.5 


-0.32 


4.6 


1.279 


1.0 


-0.32 


1.0 


















42.0-50.0 


45.78 


0.265 


0.5 


-0.32 


4.2 


1.284 


1.1 


-0.32 


1.0 


















50.0-58.5 


54.04 


0.275 


0.6 


-0.35 


3.9 


1.263 


1.1 


-0.35 


1.0 


















58.5-68.5 


63.25 


0.283 


0.6 


-0.33 


3.7 


1.261 


1.2 


-0.33 


1.0 


















68.5-81.5 


74.64 


0.283 


0.6 


-0.29 


3.6 


1.271 


1.2 


-0.29 


1.0 


















81.5-100 


90.14 


0.290 


0.6 


-0.26 


3.4 


1.281 


1.2 


-0.26 


1.0 


















100-132 


114.6 


0.290 


0.5 


-0.20 


3.3 


1.256 


1.0 


-0.20 


1.1 


















132-200 


161.3 


0.292 


0.5 


-0.22 


3.2 


1.280 


1.1 


-0.22 


1.4 


















200-300 


243.1 


0.271 


0.8 


-0.28 


3.2 


1.297 


1.6 


-0.28 


3.1 


















300-500 


382.2 


0.244 


1.2 


-0.31 


3.3 


1.428 


2.5 




7.4 




















500-1000 


687.5 


0.191 


2.2 


-0.36 


4.0 




























1000-3000 


1594 


0.087 


7.6 




8.4 











Table 6: Muon flux and charge ratio for 0.750 < cos6l < 0.825. 



momentum 


(p) 


* ■ {P? 


^|tat 


p* 




R 


A sta,t 


PR. 




interval 


[GcV] 


r GeV^ 1 
1 cm-^ s sr J 


[%] 




[%] 




[%] 




[%] 


[GeV] 




















34.5-42.0 


38.02 


0.253 


0.6 


-0.38 


4.6 


1.273 


1.2 


-0.38 


1.8 


















42.0-50.0 


45.78 


0.265 


0.6 


-0.35 


4.1 


1.292 


1.1 


-0.35 


1.8 


















50.0-58.5 


54.04 


0.275 


0.5 


-0.38 


3.7 


1.263 


1.1 


-0.38 


1.8 


















58.5-68.5 


63.25 


0.281 


0.6 


-0.34 


3.4 


1.273 


1.2 


-0.34 


1.8 


















68.5-81.5 


74.64 


0.291 


0.6 


-0.31 


3.1 


1.263 


1.2 


-0.31 


1.8 


















81.5-100 


90.15 


0.293 


0.6 


-0.27 


2.9 


1.260 


1.2 


-0.27 


1.8 


















100-132 


114.6 


0.298 


0.5 


-0.22 


2.7 


1.266 


1.0 


-0.22 


1.8 


















132-200 


161.4 


0.298 


0.5 


-0.23 


2.5 


1.281 


1.0 


-0.23 


1.9 


















200-300 


243.2 


0.285 


0.7 


-0.29 


2.4 


1.267 


1.6 


-0.29 


2.8 


















300-500 


382.3 


0.256 


1.1 


-0.31 
-0.35 


2.6 


1.394 


2.5 




5.7 




















500-1000 


688.2 


0.203 


2.0 


3.5 










1000-3000 


1601 


0.120 


5.3 




7.8 











Table 7: Muon flux and charge ratio for 0.675 < cos 6* < 0.750. 
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momentum 


<p> 


* ■ (P)^ 




p<s> 




R 


A stat 


PR 




interval 


[GeV] 


r GeV^ 1 
[ cm-^ s Sr J 


[%] 




[%] 




[%] 




[%] 


[GeV] 




















42.0-50.0 


45.79 


0.256 


0.6 


-0.41 


4.7 


1.262 


1.3 


-0.41 


1.4 






50.0-58.5 


54.05 


0.269 


0.6 


-0.42 


4.3 


1.284 


1.2 


-0.42 


1.4 


















58.5-68.5 


63.26 


0.282 


0.6 


-0.36 


4.0 


1.261 


1.3 


-0.36 


1.4 


















68.5-81.5 


74.65 


0.290 


0.6 


-0.31 


3.7 


1.277 


1.3 


-0.31 


1.4 


















81.5-100 


90.16 


0.299 


0.6 


-0.26 


3.5 


1.267 


1.2 


-0.26 


1.4 


















100-132 


114.6 


0.307 


0.5 


-0.22 


3.3 


1.287 


1.0 


-0.22 


1.4 


















132-200 


161.4 


0.312 


0.5 


-0.24 


3.1 


1.274 


1.0 


-0.23 


1.4 


















200-300 


243.2 


0.300 


0.8 


-0.32 


3.0 


1.278 


1.6 


-0.31 


2.1 


















300-500 


382.5 


0.277 


1.2 


-0.35 
-0.38 


3.1 


1.317 


2.4 




4.7 




















500-1000 


688.5 


0.231 


2.0 


3.8 










1000-3000 


1602 


0.117 


6.3 




8.7 











Table 8: Muon flux and charge ratio for 0.600 < 008 6* < 0.675. 



momentum 


<p> 


* ■ (pf 


A stat 


P4 




R 


A stat 
^R 


PR 




interval 


[GcV] 




[%] 




[%] 




[%] 




[%] 


[Gov] 




















42.0-50.0 


45.79 


0.256 


1.4 


-0.49 


5.6 


1.281 


2.8 


-0.50 


1.4 


















50.0-58.5 


54.05 


0.255 


0.9 


-0.45 


5.3 


1.269 


1.7 


-0.45 


1.4 


















58.5-68.5 


63.26 


0.273 


0.8 


-0.38 


5.1 


1.290 


1.6 


-0.37 


1.4 


















68.5-81.5 


74.65 


0.284 


0.8 


-0.32 


4.9 


1.280 


1.5 


-0.32 


1.4 


















81.5-100 


90.16 


0.294 


0.7 


-0.25 


4.7 


1.278 


1.4 


-0.25 


1.4 


















100-132 


114.6 


0.305 


0.5 


-0.22 


4.6 


1.258 


1.1 


-0.21 


1.4 


















132-200 


161.5 


0.316 


0.6 


-0.22 


4.5 


1.254 


1.2 


-0.22 


1.6 


















200-300 


243.4 


0.315 


0.8 


-0.30 


4.4 


1.285 


1.7 


-0.30 


2.6 


















300-500 


382.7 


0.297 


1.3 


-0.33 


4.4 


1.332 


2.6 




5.6 




















500-1000 


689.3 


0.243 


2.3 


-0.36 


4.8 




























1000-3000 


1604 


0.153 


5.9 




9.4 











Table 9: Muon flux and charge ratio for 0.525 < 008 6* < 0.600. 
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Figure 1: Schematic view of the experimental setup 
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Figure 2: Curvature difference at fOO GeV. The line denotes a fit with a sum 
of two Gaussian distributions with width ai and (T2. The fraction of events 
with width ai is denoted by r. 
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Figure 3: L3+C detector performance: (a) relative momentum resolution as a 
function of the muon momentum at the detector-level, (b) detector exposure 
for this analysis as a function of the muon momentum at surface for positive 
and negative muons (the sum over all zenith angle bins is shown). 



26 




200 



300 



400 500 600 700 
time since 1.1.1999 [d] 



Figure 4: Relative rate change with time for muon momenta between 50 and 
62 GcV, compared to a prediction of the atmospheric effect obtained with 
the TARGET air-shower simulation. The value of a comparison of data and 
Monte Carlo is also shown. 
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Figure 5: Relative uncertainties of the vertical zenith angle bin measurements 
for (a) the momentum scale, (b) the muon flux and (c) the charge ratio. The 
individual contributions are added in quadrature. 
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Figure 6: Momentum distribution of the selected Z-events and the back- 
ground. The Monte Carlo [31-33] events are normalized to the Standard 
Model expectation as given in Equation (11). The arrow indicates the low 
momentum cut. 
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Figure 7: The measured muon flux for zenith angles ranging from 0° (bot- 
tom) to 58° (top). The inner bars denote the statistical uncertainty, the 
full bars show the total uncertainty. For better visibility, an offset of 
0.05 GeV^cm~^s~^sr~^ was added consecutively and lines are shown to guide 
the eye. 
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Figure 8: The L3+C vertical muon spectrum compared to previous direct 
measurements providing an absolute flux normalization. All data are extra- 
polated to sea level. 
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Figure 9: The measured muon charge ratio for zenith angles ranging from 
0° (bottom) to 58° (top). The inner bars denote the statistical uncertainty, 
the full bars show the total uncertainty. For better visibility, an offset of 0.5 
was added consecutively and lines are shown to guide the eye. 
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